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Abstract 

We investigate the renormalization of the twist-3 operators which are rel- 
evant for the spin-dependent structure function g 2 (x,Q 2 ). We derive the 
anomalous dimension for the non-singlet part by calculating the off-shell 
Green's functions of the twist-3 operators including the operators which 
are proportional to the equation of motion. 



hep-ph/9603"377 



* Supported in part by the Monbusho Grant-in- Aid for Scientific Research No. 050076. 
^Supported in part by the Monbusho Grant-in- Aid for Scientific Research No. C-06640392. 



The nucleoli spin structure observed in the deep inelastic scattering is described by 
the two structure functions g\(x,Q 2 ) and g2(x,Q 2 ). In the framework of the operator 
product expansion and the renormalization group, not only the twist-2 operators but 
also the twist-3 operators contribute to g<i in the leading order of 1/Q 2 0. The general 
feature characteristic to the higher-twist operators is the occurrence of the operators 
which are proportional to the equation of motion (EOM operators) 0. Analyzing the 
twist-3 operators must be interesting and useful from the theoretical viewpoint, since 
they are the simplest examples of the higher-twist operators. 

There have been a lot of works |3], [|, || on the Q 2 -evolution of g<i- We here mention 
the two works which are frequently cited, (i) Ji and Chou || computed anomalous 
dimensions of the twist-3 operators for g 2 in the Feynman gauge (for general spin 
n). They employed the massless on-shell scheme to compute the three-point function. 
However, it is not clear how the EOM operators are dealt with since their scheme itself 
might not be consistent due to the infrared singularities coming from the collinear 
configurations, (ii) Bukhvostov, Kuraev and LipatovQ] derived GLAP-type evolution 
equation. But, this was carried out in the axial gauge and the relation to the covariant 
approach is unclear. 

The results of these two works seem not to be identical. Thus, the computation of 
the anomalous dimensions in a covariant gauge in a fully consistent scheme is desir- 
able. To do this, we decided to compute the off-shell Green's functions to renormalize 
the operators. In this case, the EOM operators should be included as independent 
operators. Infrared cut-off is provided by the external off-shell momenta. Recently this 
scheme has been employed in Ref. I| and the complete calculation of the anomalous 
dimensions for the lowest (n = 3) moment was demonstrated. The consistency and 
the efficiency of the method were also argued. Here we extend the computation to the 
case for the general moment n. We consider the flavor non-singlet case. 
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Phenomenologically, the first data of gi via the polarized deep inelastic scattering 
have been reported recently |7j]. The extensive study will be performed in HERMES. 
The theoretical determination of the Q 2 - dependence of g^ is indispensable to extract 
the physical information from experimental data. Also, the Q 2 -dependence itself will be 
checked and provide a deeper test of QCD. In view of these, it is extremely important 
to establish the theoretical prediction, which is yet controversial as discussed above. 

We follow the convention of ref . || . Let us first list up the gauge invariant twist-3 
operators which contribute to the moment /g 1 dxx n ~ 1 g2(x, Q 2 ). In the following expres- 
sions, we suppress the flavor matrices Aj for the quark field ip. 



jn-l 



Riy^- 1 = [(n - 1)V757 CT £ {W ■ • ■ 

n-l 

- tpls^D^D 1 " 1 ■ ■ ■ • ■ ■ D^-^ip] - traces , (1) 



1=1 

1 



KT'^- 1 = —m-Vn^ + Ut + U^) (Z = l,..., 71 -2) , (2) 
i^w-^-i = i n - 2 Sm^ l5l a D^ ■ ■ ■ D^Y 71 - 1 ^ ~ traces , (3) 



ryafj,i---fj, n -l ■n~2 , ~ " Q 
n n,E — 1 <-> 



Ti — 1 



2rz 

+ i)^ a D^ ■ ■ ■ £)A*»-a 7 P»-i {%jp - m^ip] - traces , (4) 



where 



Vi = i n gSifj-f 5 D' 11 ---G^ 1 ■■■D^"~ 2 Y n ~ 1 ^- traces, 
Ui = i n ~^gS^D^ ■■■G al11 ■■■D^- 2 1 ^ 1 %Ij - traces. 

In the above equations, (- traces) stands for the subtraction of the trace terms to make 
the operators traceless and is the covariant derivative. { } means symmetriza- 
tion over the Lorentz indices, S the symmetrization over /ij and g the QCD coupling 
constant, represents the quark mass (matrix). The operators in (Q) contain the 
gluon field strength G^ and the dual tensor G^ u = \e^ ua pG a ^ explicitly. The above 
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operators are not all independent but satisfy the following equation, 

n — 1 n ~ 2 

KT^- 1 = — K^- 1 + E(« - 1 - W"^- 1 + W" 1 • (5) 

n z=l 

Therefore we can exclude one operator among (|l])-(f|) to form an independent basis. 
A convenient choice of the independent operators will be (HD and It is to 
be noticed that, for the n-th moment, n gauge- invariant operators participate in the 
renormalization. 

We multiply the operators by the light-like vector A w to symmetrize the Lorentz 
indices and to eliminate the trace terms. We then embed the operators A • R a n l = 
A w ■ ■ ■ A^R^ 1 '" 1 ""- 1 into the three-point function as (0|TA-i?^(0)A M (x)^M^)|0) 
and compute the one-loop corrections. We employ the Feynman gauge and renormalize 
the operators in the MS scheme. To perform the renormalization consistently, we keep 
the quark and gluon external lines off-shell; in this case, the EOM operator mix through 
the renormalization as a nonzero operator. 

One serious problem in the calculation is the mixing of the (many) gauge non- 
invariant EOM operators. As explained in Ref . |J , these operators are given by replac- 
ing some of the covariant derivatives _D Mi by the ordinary derivatives <9 Mi in (|J). This 
problem can be overcome by introducing the vector Q^, which satisfies A M f2 M = 00, 
and by contracting the Green's functions (0|TA • l (0)A fl (x)il}(y) , 4>(z)\0) with this 
vector. This brings the two merits: Firstly, the tree vertices of the gauge invariant 
and non-invariant EOM operators coincide. Thus, essentially only one EOM opera- 
tor is now involved in the operator mixing. Secondly, the structure of the vertices 
for the twist-3 operators are simplified extremely, and the computation becomes more 
tractable. 

After the contraction with the tree vertices for A • R° nl and A • R a nE become 
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respectively, 

K,l = ((p - <7) n ~ 2 ~y- V $75 q) l - l p n - 2 - 1 $ 7 CT 75 t a , (6) 

and 

& n = ((P ~ V $75 £ + p n ~ 2 7 " 75 t a • (7) 

Here p and g are the momenta of the incoming quark and the outgoing gluon (p = 
A ■ p , q = A • q), and t a is the color matrix normalized as Tr(t a t b ) = \5 ab . 

We present the results of the one-loop calculation for the one-particle-irreducible 
three-point function with the insertion of A-R^ ll (Fig. 1). We neglect the contribution 
of the operator proportional to the quark mass, ([3]), for the moment. In the following 
expressions, Cp = (N 2 — 1)/2N C , Cq = N c , where we assume N c color, Sj = Yjk=\ 
and e = (4 — D)/2 with D the space-time dimension. 

Diagram (A) gives 

<jl -{2C F -c g )Vy: t^'^ Kk + -^-rK, 



n-2 / -i\l-k 



^i(»-i)- 2 ci n ' fc 



Diagram (B) + diagram (C) gives 



4 2 



167T 2 e 



lk=l 



+ ^2C F (2-S l -S r ^- 1 )-j^'R^ 



(n - n-fc-lCn-i (n-Z)(Z + l) 



i / c\/~i n \ "sr- ( — ■*-)' * n-i-zC^x n — 21 — 1 

htt+x {l + l)kCi+i [n - l)(n -/)(/ + 1) 

Diagram (D) + diagram (E) gives 



(9) 



167r 2 e v * 7 /(/ + 1)(Z + 2) n ' k \l(l + l)(l + 2) n-l 

n-2 I 1 \n—k (~i i ) 

+ kkl ^ ^ (n-m + m + 2) £ - ■ (10) 



Diagram (F) + diagram (G) gives 



' 2-1 



167r 2 e 



Co £ 



1 + 3 (l-2)(l-k + l) 1 



{\2l(! + l) 21(1 + !)(/ + 2) n-Z, 



7^ 



n,k 



+ 



/ + 3 



+ 



Z-2 



1 



n.l 



- £ 



21(1 + 1) 2l(l + l)(l + 2) 2(n-l) 
"-' 2 - jfe-1 _ n-2/-2 



n 



fc=m 2(n-Z-l)(n-0 
Diagram (H) gives 



7£ 



it . k 



(n-l)(n -l)(l + 2) 



(11) 



g- 



167r 2 e 



G 



' 2-1 

£ 



i 



i 



^\l-k I 21(1 + 1) 

1 



7?°" 



+ ( 1 — ^ — S'n-J-l 
1 



1 



n-2 



+ E 



2(Z + 1) 2(n-Z) 
1 n — k — 1 



/v n,f 



k — I n — Z 



(12) 



N+i \" * '* - - 2(n-Z-l)(n-Z) 
The contributions from (|3|) are easily calculated by considering the quark two-point 
Green's function ||. The renormalization constants are determined in the MS scheme. 
We summarize the final result in the following matrix form: 



/ Rn,l ' 






Zlra 






Rn,m 


L- 







I) 


Rn,m 


\ Rn,E / 














n-2), 



(13) 



where the suffix R(B) denotes renormalized (bare) quantities. We express Z^ as 



J 'j 



Si 



9 y 
16vr% ij ' 



n — 2,m, E) . 



(14) 



The relevant components of Xy are given as follows: 
U + 1)0' + 2) 



a 



G 



(l + l)(l + 2)(l-j) 

i +j n- 2 Q-i (n-l + l- ]) 



+ (2C F - C G ) 



■1) 



2(-iy 



,_ 2 CU (n-l)(Z-j) 
(1<J<Z-1), 
1 1 



Z(Z + l)(Z + 2) 



iCj 



(15) 



x. 



a 



G~ 



+(2C F - C G ) 



(n - l-j)(n-j) 



1 



+ 



2(-l) 1 



n-l Z(Z + l)(Z + 2) n-l 
+Cf (3 — 25; — 2S n ^i-\) , 



(16) 



n-l-l){n-l){j-l) 

, n -2Cj(n- 1-1+ j] 



+ (2C F - C G ) 



Irn 



AC F 
nl(l + !)(/ + 2) 



v ' n . 2 Q (n-l)(j-l) 
(I + 1 < j < n - 2) , 
2C F 



X, 



IE 



(n -1)(Z + 1)(Z + 2) 



+ (-1 



X„ 



^ n -j n-2-lC n -2-j 
fl — I 



-AC F S,, 



n-l 



(17) 
(18) 



For the physical quantity (moments), only the Z\j , Zi m , Z mm components of the renor- 
malization matrix give the contributions since the physical matrix element of the EOM 
operators turns out to be zero. With the above X^, the anomalous dimension matrix 
for the twist-3 operators becomes, 



7ij 



^X 



(19) 



This matrix enters into the renormalization group equation for the Wilson's coefficient 
function associated with the corresponding operator as, 





9 a^ d 



E,. = 



In the present study, we have obtained the anomalous dimension for the twist-3 
operators which contribute to g 2 {x, Q 2 ). We performed the calculation with the mani- 
fest Lorentz covariance being kept. We adopted the Feynman gauge and dimensional 
regularization. We have chosen the operators which include the gluon field strength 
explicitly as an independent operator's basis although this choice of basis is never 
compulsory. To identify the renormalization constants correctly, the off-shell Green's 
functions are considered. We have shown that the EOM operators play an important 
role to complete the renormalization of composite operators and the structure of the 
renormalization constant matrix takes the triangular form expected from the general 
argument [|]. 



If we could calculate the "on-shell" matrix elements of composite operators in terms 
of purely perturbative Feynman graphs, we could obtain the enough informations with- 
out considering the EOM operators. However the infrared singularity coming from the 
collinear configuration can not be regulated ||10|| . Therefore we believe that calculating 
the off-shell Green's functions is the safest method to obtain the anomalous dimensions. 

As a technical issue, we have used the projection introduced in Ref.|| to avoid the 
complexity stemmed from the fact that there are a lot of EOM operators which are 
not gauge invariant. This projection may have some relation to the case in which the 
light-cone gauge fixing A^A^ = is adopted since this gauge does not discriminate 
between the covariant and the partial derivative in the composite operators. 

Since the authors of Ref. || adopted the on-shell scheme to extract the anomalous 
dimensions, we can not compare our results with theirs graph by graph. However 
our final results are not in agreement with theirs although our definition of composite 
operators and choice of gauge are the same as those in Ref. [0] . If one replaces the index 
i by n — 1 — i in their results, the both results become the same. 

On the other hand, our final results agree with those of Ref. @]. The approach 
adopted by them is quite different from ours. The main differences are : (a) they con- 
sidered the physical quantities (structure functions) from the beginning, so the EOM 
operators do not appear explicitly, (b) they adopted the axial gauge to derive the 
GLAP-type evolution equation. In their analyses, the evolution of the quark bilinear 
operator ([!]) was also considered. Therefore it seems that their choice of the operator 
basis is different from ours at the intermediate stage of calculations. By including the 
quark bilinear operator ([I]), we can eliminate the gauge invariant EOM operator from 
the independent operator basis (see (|)). Furthermore, they included some contribu- 
tions from the "one-particle-reducible" diagrams to take into account the circumstance 
that the partons can not be regarded as located on the mass shell |3j]. This may cor- 
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respond to the fact that even if the gauge invariant EOM operator is discarded, the 
gauge non- invariant EOM operators turn out to mix through the renormalization ||. 

We expect that future precise measurements on g 2 will clarify the effect of twist-3 
operators which may be the first quantity to see the higher- twist effect in QCD. 



We would like to thank H. Kawamura and Y. Koike for valuable discussions. 
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Figure Captions 
Fig. 1 One-loop corrections to A • R a nl . 
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